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Abstract 

Piecewise  linear  maps  over  triangulations  are  used  extensively  in  geometric  modelling  and 
computer  graphics.  This  short  note  surveys  recent  progress  on  the  important  question  of 
when  such  maps  are  one-to-one,  central  to  which  are  convex  combination  maps. 


1 Introduction 

Piecewise  linear  maps  over  triangulations  have  several  applications  in  geometric  model- 
ling and  computer  graphics.  For  example,  Figure  la  shows  a surface  triangulation  T of  a 
set  of  points  ( Xi,yi , Z{ ) sampled  from  some  unknown  surface  in  IR'\  A standard  approach 
to  fitting  a smooth  parametric  surface  s(u,  v ) to  these  points  is  to  first  parameterize  them, 
i.e.,  compute  planar  points  ( ) corresponding  to  the  data  points  (xi,yi,Zi).  Then  us- 
ing some  scattered  data  method,  we  find  a parametric  surface  s : fi  — ► It  ',  defined  over 
some  suitable  domain  0 containing  the  points  (uj,v ,),  such  that 

s(uitVi)  « (Xi,pi,Zi). 

A choice  of  parameterization  is  shown  in  Figure  lb  and  a least  squares  surface  approx- 
imation using  bicubic  B-splines  is  shown  in  Figure  lc. 

Notice  that  the  choice  of  parameter  points  ( Ui,Vj ) uniquely  determines  a piecewise 

linear  map  <j>  : Dq > IR2,  where  Dr  is  the  union  of  the  triangles  in  T.  In  practice, 

a necessary  requirement  on  0 to  ensure  adequate  quality  of  the  subsequent  surface  ap- 
proximation s(u,  v)  is  that  <j>  should  be  injective.  In  Figure  lb  the  mapping  0 was  taken 
to  be  a so-called  convex  combination  map,  which,  as  we  will  see  later,  is  guaranteed  to 
be  one-to-one  since  the  boundary  of  T is  mapped  to  a rectangle.  Put  another  way,  none 
of  the  triangles  in  Figure  lb  are  ‘folded  over’.  In  fact  further  properties  of  the  map  are 
important,  such  as  linear  precision,  and  this  was  achieved  in  Figure  lb  by  using  the 
so-called  shape-preserving  weights  (the  coefficients  in  the  convex  combinations).  For  a 
discussion  of  that,  see  [3]. 

Another  application  of  piecewise  linear  maps  is  to  image  morphing.  Image  morphing 
can  be  carried  out  by  continuously  transforming  one  planar  triangulation  T°  (whose 
vertices  represent  feature  points  in  the  image)  to  another,  T1.  Here  we  assume  that 
there  is  a one-to-one  correspondence  between  the  vertices,  edges,  and  triangles  of  T°  and 
T1.  We  can  view  each  intermediate  triangulation  T (f),  0 < t < 1,  (where  T(0)  = T° 
and  T(l)  = T1)  as  the  image  of  a piecewise  linear  map  < j>{t)  : Dr<>  — > Dr(t )■  As  with 
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parameterizations,  it  is  again  important  that  0(f)  is  one-to-one.  Figure  2 shows  a so- 
called  convex  combination  morph  of  [4]  of  two  given  planar  triangulations:  T°  appears 
on  the  left  and  T1  on  the  right.  The  two  triangulations  in  the  middle  are  T{  1/3)  and 
T( 2/3).  This  morph  ensures  that  (f>{t)  is  one-to-one  for  all  t in  [0, 1]  and  therefore  T(t) 
has  no  ‘folded’  triangles  at  any  time  instant  t. 

Piecewise  linear  maps  also  arise  in:  texture  mapping;  numerical  grid  generation;  and 
in  setting  up  multiresolution  frameworks  (nested  spaces  of  piecewise  linear  functions) 
for  manifold  surface  triangulations  in  computer  graphics. 


Fig.  1.  Spatial  triangulation  (la),  Convex  combination  parameterization  (lb),  Bicubic 
spline  approximation  (lc). 


Fig.  2.  Convex  combination  morph. 


2 Convex  combination  maps 

For  the  sake  of  simplicity  we  will  only  discuss  convex  combination  maps  defined  over 
planar  triangulations  even  though  all  the  results  hold  equally  well  when  the  domain  of  the 
map  is  a spatial  triangulation  such  as  that  in  Figure  la.  Thus  let  T = {TV, . . . , Tm}  be  a 
simply-connected  planar  triangulation,  with  closed  triangles  7),  and  let  Dr  = Utgt^1’ 
as  in  Figure  3.  We  will  call  a mapping  0 : Dq — > R2  a convex  combination  map  if 
it  is  piecewise  linear  over  T and,  for  every  interior  vertex  v of  T,  there  exist  weights 
A vw  > 0,  for  w £ Nv,  such  that 

^ ' AyW  = 1, 
weN„ 
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and 

(p{v)  = 5Z  (1) 

t«€Ar„ 

where  Nv  is  the  set  of  neighbours  of  v;  see  Figure  3. 


Fig.  3.  Convex  combination  map. 


In  applications,  the  mapped  boundary  vertices  f>{v)  are  chosen  first.  Then  the  weights 
\vw  are  all  specified  according  to  some  chosen  strategy.  Then  finally  the  mapped  interior 
vertices  are  found  by  treating  the  equations  in  (1)  as  a linear  system. 

Example  2.1  If  an  interior  vertex  v of  T has  five  neighbours  then  we  might 

set  11111 

<f>(v)  = + ~4>{v2)  + gt(v3)  + j0(v4)  + 

Until  recently,  the  only  theory  behind  convex  combination  maps  was  that  of  Tutte  [8]. 
Working  from  a purely  graph-theoretic  point  of  view,  Tutte  proposed  a so-called  bary- 
centric mapping  for  constructing  straight  line  drawings  of  3-connected  graphs  (which 
include  triangulations).  A barycentric  mapping  in  our  context  is  simply  a convex  com- 
bination map  in  which  all  the  weights  at  each  vertex  are  equal,  i.e.,  Xvw  = l/dv,  where 
dv  is  the  degree  or  valency  of  the  vertex  v.  Thus  for  v in  Example  1 we  must  have 

l)  + g<A(v2)  + + g0(u4)  + ^4>(v  5). 

Tutte  showed  that  a valid  straight  line  drawing,  i.e.  one  with  no  edge  crossings,  results 
from  a barycentric  mapping  if  the  ‘boundary’  of  the  graph,  a so-called  ‘cycle’,  is  mapped 
to  a convex  polygon.  However,  as  argued  in  [3],  convex  combination  maps  share  all  those 
properties  of  barycentric  maps  necessary  for  Tutte’s  proof.  Thus  when  interpreted  in  the 
right  way  and  suitably  generalized,  Tutte’s  theorem  can  be  expressed  in  the  following 
way. 

Theorem  2.2  Suppose  <j> : Dq — > R2  is  a convex  combination  mapping  which  maps  the 
n boundary  vertices  ofT  cyclically  into  the  n vertices  of  some  n- sided  convex  polygon  in 
the  plane.  Then  cf>  is  one-to-one. 

Despite  this  generalization,  however,  there  are  still  two  aspects  of  it  which  need  to 
be  improved  from  the  point  of  view  of  applications  and  future  research. 

The  first  is  that  we  would  like  to  extend  the  theorem  so  that  we  can  allow  some,  and 
indeed  many,  of  the  mapped  boundary  vertices  to  be  collinear.  Indeed  in  the  application 
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to  parameterization  for  surface  fitting,  it  might  be  convenient  to  map  all  the  boundary 
vertices  of  the  given  triangulation  into  the  four  sides  of  a rectangle,  as  in  Figure  lb. 
This  is  because  tensor-product  splines  surfaces  are  defined  over  rectangular  domains. 
Collinearity  will  also  often  be  desirable  in  morphing,  as  in  Figure  2,  and  in  most  other 
applications.  Thus  a drawback  of  Theorem  2.2  is  that  it  does  not  allow  collinear  vertices 
in  the  image  boundary. 

The  second  aspect  is  that  we  would  like  to  simplify  the  proof  in  order  to  have  some 
hope  of  establishing  the  injectivity  of  piecewise  linear  maps  in  even  more  general  situ- 
ations, such  as  when  mapping  to  non-convex  regions,  or  when  some  of  the  mapped  ver- 
tices are  constrained,  for  example.  The  fact  that  Tutte’s  proof  relies  on  the  non-existence 
of  the  Kuratowski  subgraphs  K$  and  3 in  a planar  graph  illustrates  its  complexity. 

It  is  these  two  improvements  that  are  the  focus  of  [5].  The  main  idea  of  [5]  is  the 
observation  that  Theorem  2.2  is  very  similar  to  a theorem  on  harmonic  maps,  referred 
to  by  Duren  and  Hengartner  [2]  as  the  Rado-Kneser-Choquet  thereom,  which  was  es- 
tablished in  [7,  6,  1].  Recall  that  a mapping  <f> : D — ► 1R2,  with  D C 1R2  and  <f>  = (u,  v), 
is  harmonic  if  both  its  components  u(x,y)  and  v(x,y)  satisfy  the  Laplace  equation  in 
D,  i.e., 

'U’yy  ~~  d,  VXx  T V yy  — Oj 

see  Figure  4. 

Rado-Kneser-Choquet  Theorem.  Suppose  </>:£>—>  R2  is  a harmonic  mapping 
which  maps  the  boundary  dD  homeomorphically  into  the  boundary  8Q  of  some  convex 
region  QcR2.  Then  <j>  is  one-to-one. 


Fig.  4.  Harmonic  map. 


This  suggested  that  a proof  of  Theorem  2.2  might  be  based  on  a proof  of  the  Rado- 
Kneser-Choquet  theorem,  in  particular  the  short  proof  of  Kneser  [6].  Kneser’s  proof 
begins  by  showing  that  <j>  is  locally  one-to-one  in  the  sense  that  the  Jacobian  of  (p, 

Ux  'U>y 
VX  Vy 

never  vanishes.  Kneser  establishes  this  by  supposing  that  the  Jacobian  is  zero  at  some 
point  (a;0,  yo)-  In  that  case  there  must  be  a straight  line  ax + by  + c = 0 passing  through 
the  point  <f>(xo,  yo)  such  that  both  partial  derivatives  of  the  function  h(x,  y)  = au(x,  y)  + 
bv(x,y)  + c are  zero  at  (xq, yo).  At  the  same  time,  the  function  h : D — > R is  zero  at 
(xq,  yo)  and  has  just  two  zeros  along  the  boundary  of  D.  Noting  that  h(x,  y)  is  a harmonic 
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function,  Kneser  then  uses  the  Nodal  Lines  theorem  of  Courant  to  argue  that  there  are  at 
least  four  zero  contours  of  h emanating  from  (x0.  y0)  and  due  to  the  maximum  principle 
for  h,  these  four  curves  can  never  self-intersect,  nor  intersect  one  another.  Therefore  all 
four  curves  must  reach  the  boundary  of  D which  is  a contradiction. 

These  ideas  were  used  in  [5]  to  establish  a much  simpler  proof  of  Theorem  2.2  than 
that  of  Tutte.  No  graph  theory  is  needed  at  all.  Instead,  the  discrete  maximum  principle 
for  convex  combination  functions  plays  the  role  of  the  maximum  principle  for  harmonic 
functions.  Similar  to  Kneser’s  proof  we  show  first  that  </>  is  locally  one-to-one,  except 
that  we  understand  this  to  mean  that  the  restriction  of  (j>  to  any  quadrilateral  in  T is 
one-to-one,  a quadrilateral  being  the  union  of  two  triangles  sharing  a common  edge. 


Vl 


Moreover,  Theorem  2.2  is  generalized  in  [5]  to  allow  collinear  mapped  boundary 
vertices.  We  call  an  edge  [u,  w]  of  T a dividing  edge  if  both  endpoints  v and  w are 
boundary  vertices  yet  the  edge  [«,  w]  itself  is  not  contained  in  the  boundary.  For  example 
in  Figure  5,  the  only  dividing  edge  in  the  triangulation  is  [tq,^].  Dividing  edges  play  a 
critical  role  because  they  partition  the  triangulation  into  subtriangulations  %,  in  each 
of  which  every  convex  combination  function  satisfies  a discrete  maximum  principle  in  its 
strong  form.  The  main  result  of  [5]  was  the  following. 

Theorem  2.3  Suppose  T is  any  triangulation  and  that  <f>  : Dj  — > R2  is  a convex 
combination  mapping  which  maps  dD?  homeomorphically  into  the  boundary  dQ  of  some 
convex  region  0 C R2.  Then  <p  is  one-to-one  if  and  only  if  no  dividing  edge  \v,w\  of  T 
is  mapped  by  cp  into  dfl. 

3 Future  research 

Here  is  a list  of  topics  for  future  research. 

• A triangulation  is  a special  (maximal)  kind  of  planar  graph.  Can  one  extend  The- 
orem 2.3  to  other  planar  graphs,  for  example,  rectangular  grids?  This  is  likely 
because  Tutte’s  theory  already  holds  for  all  3-connected  graphs. 

• In  what  way  can  the  theorem  be  extended  from  bivariate  maps  to  trivariate  ones? 

• Can  similar  one-to-one  maps  be  guaranteed  when  mapping  closed  surfaces  of  vari- 
ous topology?  For  example,  we  would  like  to  map  a closed  manifold  triangulation, 
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homeomorphic  to  a sphere,  into  a unit  sphere  injectively.  Here  each  triangle  in  the 
triangulation  would  be  mapped  to  a spherical  triangle  on  the  surface  of  the  sphere. 

• Can  one  find  sufficient  conditions  for  the  injectivity  of  constrained  maps,  i.e. , piece- 
wise  linear  maps  in  which  the  image  of  certain  interior  points  is  specified  in  advance? 

• Can  one  remove  the  requirement  of  having  to  map  the  boundary  to  a convex  polygon 
and  still  ensure  a one-to-one  mapping  under  some  weaker  condition? 

• Can  the  Rado-Kneser-Choquet  thereom  and  Theorem  2.3  be  combined  as  part  of 
a single  more  general  theorem? 
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